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The Mandelstam-Tamm and Margolus-Levitin inequalities play an important role in the study of quantum me- 
chanical processes in Nature, since they provide general limits on the speed of dynamical evolution. However, to 
date there has been only one derivation of the Margolus-Levitin inequality. In this paper, alternative geometric 
derivations for both inequalities are obtained from the statistical distance between quantum states. The inequal- 
ities are shown to hold for unitary evolution of pure and mixed states, and a counterexample to the inequalities 
is given for evolution described by completely positive trace-preserving maps. The counterexample show that 
there is no quantum speed limit for non-unitary evolution. 



I. INTRODUCTION 

In experiments, we often do not have direct access to the pa- 
rameter we want to measure, but only to certain observable 
variables. We tend to formulate a model of the system that 
tells us how we expect the observable variables depend on 
the parameters of interest of the system. Given an estimator 
function based on these observables, we can then calculate the 
most likely value of the parameter based on the measurement 
data [1-4]. When estimating an unknown parameter we will 
always be restricted by statistical uncertainty: the more tests 
we perform the greater our confidence in the result. More 
quantitatively, given a large number of independent tests A^, 
the average error is at best AO = a /^/N, where is the vari- 
ance in each test. In quantum metrology it is possible to beat 
this limit by using entanglement. The fundamental Heisen- 
berg limit then allows us to estimate the parameter Q with 
average error A0 °^ 1 /N, where the tests are no longer inde- 
pendent of each other In order to accommodate entanglement 
across the tests, the estimation procedure must essentially 
boil down to a single-shot measurement. To distinguish be- 
tween two configurations in a single-shot measurement with 
unit probability, the two configurations must be described by 
orthogonal quantum states [5-7]. If Q parametrizes the evolu- 
tion of the initial state to the final state, maximum sensitivity 
is obtained when the state of the system evolves to an orthog- 
onal state for the smallest possible value of 9. This is the 
geometric interpretation of parameter estimation [8-1 1]. 

The minimum value of Q needed for a quantum system to 
naturally evolve to an orthogonal state is determined by the 
moments of the generator K of rotations (or translations) in 
Q. In particular there are two bounds on this evolution. First, 
for the specific case of time evolution (0 = f ) it was shown by 
Mandelstam and Tamm that a lower bound could be defined 
in terms of the variance in the energy of the system [12]: 

t> , (1) 

with (A£')2 = - {hY, and H is the Hamiltonian of the 
system. Relatively recently, Margolus and Levitin derived a 
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second bound on the the minimum time of unitary evolution, 
which is instead determined by the average energy of the sys- 
tem [13]: 



where E = {H) . For the purposes of determining the so-called 
maximum speed of evolution for a quantum system, having a 
bound defined in terms of the average energy of the system is 
in general far more useful than a bound in terms of the vari- 
ance, since the average energy is generally easier to determine 
than the variance in the energy. For example, the Margolus- 
Levitin bound was used to estimate the ultimate limits to com- 
putation and the maximum computational power of the uni- 
verse [14, 15]. 

In this paper, we give an alternative derivation of the 
Mandelstam-Tamm and Margolus-Levitin inequalities based 
on the distance of quantum states in Hilbert space. In section 
II we give a brief review of the concept of statistical distance, 
both for classical probability distributions and for quantum 
states, and we relate the statistical distance to the Fisher in- 
formation, which measures the amount of information about a 
parameter obtained in a given measurement. In section III we 
define the speed of dynamical evolution as the rate of change 
of the statistical distance, analogous to the velocity in real 
space, and we use this to derive the two inequalities. Finally, 
we will show by means of a counterexample that the inequali- 
ties apply to unitary evolutions, and can be violated in general 
completely positive trace preserving maps (CP maps). 



II. STATISTICAL DISTANCE 

In order to quantify the difference between two probability 
distributions, we can define a distance measure in the space 
of probability distributions, called the statistical distance. In 
general, a distance s between two points a and i> in a metric 
space has the following properties: 

1. s{a,b) > 0, 

2. s{a,b) =0 4^ a = b, 

3. s{a,b) — s{b,a), 

4. s{a,c) < s{a,b) +s{b,c) (the triangle inequality). 
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When the two points a and b are very close together (b = 
a + da), we can define the infinitesimal distance between the 
points in terms of a metric tensor gji;, according to 

ds^ = ^gjiida'da'^ , (3) 

where the da^ are the components of the tangent vector to 
a. We make the distinction between upper and lower indices, 
with the upper indicating the tangent (contravariant) vector, 
and the lower index indicating the one-form (covariant vec- 
tor). In the case of Euclidean spaces the metric will be pro- 
portional to the Kronecker delta Sjk, and Eq. (3) reduces to 
Pythagoras' theorem. 

A. Statistical distance for classical probability distributions 

In general, one can define many different distance measures 
on a metric space. As a trivial example, we can multiply any 
distance function by an arbitrary constant and still satisfy the 
four properties above. This corresponds to changing the units 
of length for the distance function. In our case, we want to 
choose a natural distance function for the space of probability 
distributions that relates the statistical distance to expectation 
values. The natural inner product in the space of probability 
distributions is the expectation value of a random variable A 

(A)=£A,y, (4) 
j 

and the correlation is 

(AB) = J^AjBjpJ = Y^AjB^gj' . (5) 

j jk 

We related the correlation to the statistical distance by writing 
(AB) in terms of a metric g^'^. Since the probability distribu- 
tions form the vectors in the space, the probabilities are the 
proper tangent vectors, while the one-forms A, and Bj^ are the 
values of the random variables A and B. It is clear that the 
correlation forms the natural quadratic form for the space of 
probability distributions. 

From Eq. (5), we find that the metric g^'^ must be propor- 
tional to p^, together with a Kronecker delta 5jk to match the 
indices. However, the metric we require for the infinitesimal 
form of the statistical distance in Eq. (3) has lower indices. 
We therefore find that 

gjk = {g^')-' = ^f. (6) 

Consequently, the distance between two probabihty distribu- 
tions p^ and pi + dpi using this metric becomes 

^^' = I^,xW = I^. (7) 

hk j P 

This is the infinitesimal statistical distance for classical prob- 
ability distributions. The distance increases without bound 



when one of the ph becomes zero. This is interpreted as 
follows: when we try to distinguish between two probabil- 
ity distributions and pg, and one type of event has zero 
probability in p^ but not in ps, then the occurrence of that 
event immediately tells us with certainty that our system is 
described by pB- In practical numerical applications, we ex- 
clude the boundary of the probability simplex to avoid these 
divergency issues. 

At this point, we note a peculiarity in the form of the sta- 
tistical distance for classical probability distributions. If we 
write the probability pi = {ri)-^ in terms of (real) probabil- 
ity amplitudes r ', we obtain dpi = Iridri , and the statistical 
distance becomes 

ds^=AY,{dr'f. (8) 
j 

In other words, the statistical distance in the space of probabil- 
ity distributions is the Euclidean distance in the space of (real) 
probability amplitudes. However, probability amplitudes are 
usually associated with quantum mechanics, and our entire 
discussion has been classical. Note also the factor 4 in the sta- 
tistical distance. We could remove it by rescaling the distance 
units. However, our units are quite convenient, and the factor 
will become important in section III, when we calculate the 
speed of dynamical evolution. 

So far, we have been dihgent in observing the difference 
between upper and lower indices in order to derive the correct 
form of the metric tensor. In the remainder of this paper, this is 
no longer necessary, and we will use only lower indices from 
now on. 



B. Statistical distance for density matrices 

Next, we will derive the statistical distance for density matri- 
ces p. We will follow the general procedure of the previous 
section, and derive the equivalent of the metric tensor by iden- 
tifying the natural quadratic form. Again, the natural inner 
product on the space of density matrices is the Born rule for 
the expectation value of a quantum mechanical observable A, 
given by 

(A)=Tr(pA). (9) 

Similarly, we can choose the correlation as the natural 
quadratic form. However, there is a subtlety: since corre- 
lations are observable, they must be represented by Hermi- 
tian operators. On the other hand, when we consider the cor- 
relation between two non-commuting observables the opera- 
tor product AB is not Hermitian: {AB)^ = B^A^ =BA^ AB. 
The correct correlation is therefore given by the symmetrized 
product of A and B: 

^(AB+BA) = ^Tr[p{A,B}], (10) 

where we used the anti-commutator {A,B} = AB-\-BA. Using 
the cyclic property of the trace, this can be rewritten as 

^-{AB+BA)=Tx[AMp{B)\, (11) 
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with 



= \{p,B}^^-Y,{pj+pk)Bj,\j){k\, 



and its inverse 



-'(B) 



2B 



\j){k\- 



(12) 



(13) 



Here, we used the diagonal basis for p = Jlj Pj\j) {J\- The 
"raising" and "lowering" operators and play the role 
of the metric, and the infinitesimal statistical distance between 
p and p + dp on the space of density operators can be written 
as [9] 



dsl=Tr[dp^p{dp)] 



(14) 



In the case of pure quantum states this distance simplifies con- 
siderably, and we obtain the Wootters distance [8]. If I'/) and 
are two arbitrary pure quantum states in Hilbert space, 
then the statistical distance is given by the angle between the 
two states: 



s(VA,(?)) = ai-ccos|(v/|0)|. 



(15) 



Notice how the Wootters distance is related to the fidelity 
I (1/10) which is the probability of mistaking jy/) for |0) 
in a single-shot measurement. 



C. Relation to the Fisher Information 

The statistical distance is a measure of how close one proba- 
bility distance is to another. In other words, it can be directly 
related to the number of times we need to sample our system 
in order to tell reliably which of the two probability distribu- 
tions describes our system. Imagine that the two probability 
distributions are connected by a curve that is parametrized by 
a real number 9. Distinguishing between the two distribu- 
tions then reduces to the estimation of the parameter 9. This 
is a well-known problem, and the information about in a 
particular measurement procedure is given by the Fisher in- 
formation F{9). For a discrete set of possibilities, the Fisher 
information is 



Moreover, using the expression in Eq. (14) for the quantum 
mechanical statistical distance, we obtain 



nO)={^) =Tr [p%(p')] 



(18) 



where p' is the derivative of p with respect to 0. If translations 
in 9 are generated by a Hermitian operator K, we can use the 
Heisenberg equation of motion to write 

where (K) = Tr(pA') is a real number, and can therefore be 
included in the commutator with impunity. Substituting this 
into Eq. (18) yields [16] 



F{9) ^ Tr[p'^pip')]=l^j:^-P^-^ 



(20) 



In other words, the amount of information about 9 in any mea- 
surement is bounded by the variance in its generator [16]. 



III. DYNAMICAL EVOLUTION OF QUANTUM STATES 

The quantum speed limit for the dynamical evolution of 
quantum states will take the form of two inequalities, the 
Mandelstam-Tamm inequality and the Margolus-Levitin in- 
equality. Before we derive these inequalities, we will have 
to define precisely what we mean by the speed of dynami- 
cal evolution. Analogous to the instantaneous velocity of a 
particle in real space, which is the time derivative of the posi- 
tion function of the particle, we define the dynamical speed of 
evolution for quantum states as the derivative of the statistical 
distance function with respect to the parameter 9: 



ds 
d9 



(21) 



This allows us to use the results from the previous section to 
obtain bounds on v(0). 



F(9) 



1 



'pjW 



d9 



ds 
dO 



where Pj{9) is determined by the Born rule 

Pj{9)^lr[Ejp{9)], 



(16) 



(17) 



with Ej the POVM associated with measurement outcome /'. 
A similar expression holds when the possible events form a 
continuum. 

We see from Eq. (16) that the Fisher information is directly 
related to the derivative of the statistical distance, as expected. 



A. Mandelstam-Tamm inequality 

The Mandelstam-Tamm inequality is almost immediate from 
Eq. (20): 



Taking the positive roots, we can rewrite this as 



ds 2 , 

d-9^-n'^^ 



(23) 
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where 5K = \ ^{{AKY) \ . Separating the variables and inte- 
grating yields 



L 



^ \ ti 

dd'> - ^ ds 
- 2 5kJq 



In the case where 9 is the time parameter generated by 
the Hamiltonian H, the inequality reduces to the famous 
Mandelstam-Tamm inequaUty 



n h 

t> . 

- 2AE 



(25) 



Note that the integral over ds in Eq. (24) runs from to 7t, 
instead of n/2. This is due to the factor 4 in Eq. (8). Alter- 
natively, this can be seen in the Bloch sphere, were the angle 
between orthogonal qubit states is n, rather than n/2. Note 
that the Mandelstam-Tamm inequality in Eq. (25) was derived 
for unitary evolution of arbitrary mixed states. We will see 
in Section IV that this bound can be violated by non-unitary 
evolutions of density operators. 



B. Margolus-Levitin inequality 

Instead of using the bound on the Fisher information directly 
to derive the Mandelstam-Tamm inequality, we can use the ex- 
pression for the statistical distance and carry out the differen- 
tiation with respect to 6 exphcitly. Remarkably, this will yield 
the Margolus-Levitin inequality. We will proceed by first de- 
riving the inequality for pure states, followed by an extension 
to mixed states via a standard purification procedure. 

We again consider the evolution parameterized by Q, which 
is generated by the Hermitian operator K. In this case, a sys- 
tem initially described by the pure state |v/o) at = will 
evolve to 



IV^e) =exp( --/:0 ) IvAo). 



(26) 



The Wootters distance between two pure states is given by the 
angle between the states and the rate of change of the statisti- 
cal distance can therefore be written as 



ds 
dd 



d_ 
dO 



arccos(|(v/b|V^e)|) 



^/l-\{Wo\We)\^de 



\{Wo\We)\- (27) 



Since the pre-factor 1 / ^/{^^^^ > 1 for all real x, we obtain 
the inequality 



ds d , , , ., 



Next, we prove that 

_d_ 
dd 



\{¥o\Ve)\ < 



dd 



(28) 



(29) 



To this end, we rewrite the derivative on the left-hand side of 
Eq. (29) as 

^KV^olV^e)! = ^V(V^o|V^e)(V^elV^o), (30) 
and using the generalised Schrodinger equation 

in-^WB)=K\^ife), (31) 



this becomes 
d 

dd 



\{¥o\¥e}\ = 



< 



-i{Wo\K\\i/e){Ve\¥o) + i{Vo\Ve){Ve\K\Wo) 

2h\{xi/o\\i/e}\ 
lm{{\i/o\K\\i/g){\i/e\Wo)) 

^\{¥oWe)\ 
|(yo|g|vA9)(v/e|v/o)| 

^KV^olV^e)! 



(32) 



The right-hand side of Eq. (29) becomes 
d 



dd 



= -n\<^o\K\We}\ = ^^^^^ 



> 



\{Wo\K\Ye){Vo\Wd)\ 



(33) 



where we again used the generalized Schrodinger equation, 
and in the last line we used the Cauchy-Schwarz inequality. 
Finally, we combine Eq. (33) and Eq. (32) to obtain 



dd 



|,„ M ^ M{Vo\K\Ve){Ve\Vo}) . 



dd 



(34) 

Since arccos(x) is a monotonically decreasing function in the 

interval < x < 1 , the derivative of | ( i/Zq | \/Ag ) ) | with respect to 
Q is strictly positive and we have therefore proved Eq. (29). 

Continuing the derivation of the Margolus-Levitin inequal- 
ity, we use Eq. (34) in Eq. (28) and find that 



ds 
dd 



< 



dd 



iWolWe) 



^ \{Wo\K\We)\ ^ \{Vo\K\Vo)\ 
- h - h 



Separating the variables s and d, we obtain 

h 



m\ 

and integrating both sides gives 



ds , 



2 \{K)\ 



m 

h ■ 

(35) 



(36) 



(37) 



In the case where K is the Hamiltonian and the time, the 
inequality becomes the Margolus-Levitin inequahty 



t > 



n h 
2E 



(38) 



with E the average energy of the system. 



5 



So far, we have shown that the Margolus-Levitin inequal- 
ity holds for the unitary evolution of pure states. However, 
following the derivation of the Mandelstam-Tamm inequality, 
we would like to extend this bound to the unitary evolution of 
arbitrary mixed states. The challenge is to find the most con- 
venient distance measure between mixed states, i.e., a gener- 
alization of the Wootters distance for density operators. We 
choose the purification of the density matrices that maximizes 
the fidelity between them. The Margolus-Levitin inequality 
then applies to the purifications, which in turn can be trans- 
lated to a bound on the speed of unitary evolution of mixed 
states. 

The fidelity between two density matrices p and a is de- 
fined as [17] 



Tr 



(39) 



which can be interpreted as the probability of mistaking p 
for a in a single-shot measurement. Despite its appearance, 
F{p,a) is symmetric in p and a. We can relate the fidelity to 
a distance measure in various ways, one possibility being 



s(p, cj) = arccos \/F{p,a) 



(40) 



When both p and <7 are pure states s{p,o) reduces to the 
Wootters distance of Eq. (15). According to Uhlmann's theo- 
rem [18], the square root of the fidelity is given by 



where {Ks) is the expectation value of K on the system S, 
and {Kr) the expectation value on R. In order to interpret 
this inequality, let K be the Hamiltonian of the system. The 
expectation value {Ks) + {Kr) is then the average energy of 
the system S and the purification system R taken together At 
first glance this seems problematic since the purification is 
merely a mathematical construction and therefore unrestricted 
amounts of energy may be added to the system during the pro- 
cess. This would lead to arbitrarily short orthogonality times. 
However it is precisely because the purification is a mathe- 
matical construct that this is not a problem: any extra energy 
added is physically meaningless and therefore the true bound 
on the evolution occurs for {H) = {Hs) . In other words we 
simply have to choose purifications made up only of degener- 
ate ground eigenstates (with E = 0). This means that the evo- 
lution of mixed states is also bounded by Eq. (38). A similar 
line of reasoning holds for general Hermitian operators K. It 
was show in Ref. [19] that mixed states can never attain either 
equality. In addition, any state that does attain the inequalities 
has the form [19, 20] 



\W) 



V2 



{\xifo)+e''\Wn)), 



(44) 



where the states | xj/n) are (possibly degenerate) energy eigen- 
states of the system. 



IV. NON-UNITARY EVOLUTION 



\/f{P^=. max KxplXa) 

\Xp)'\X<r) 



(41) 



where \Xp) and \Xa) are purifications of p and a. In other 
words, if p is the state of a system S, and \Xp) is a pure state 
of a compound system S + R, then p — TirHXp) {Xp\) ■ The 
system R must be described on a Hilbert space that is at least 
as large as the Hilbert space of system S. A similar definition 
holds for \Xa)- For our purposes it is sufficient to note that 
there exist states \Xp) and \Xa) for which the equality holds: 

Fip,a)^\{Xp\Xa)\^. 

For pure states the Wootters distance is interpreted as the 
angle between the states in Hilbert space. For mixed states, 
this interpretation makes sense only when the states have the 
same purity Tr(p^) = Tr(ff^), and the evolution is unitary. 
The speed of evolution can then be bounded using 



ds 
dO 



dd 



= ^ai-ccos|(xo|Ze> 



arccos ^/F{po,Pe) 

uu 

(42) 

where po is the initial state, and pe is the evolved state. The 
pure states \xo) and \xe) are the respective purifications that 
maximize the fidelity. Since \xo) and \xe) are pure, this is 
just the derivative of the standard Wootters distance in the 
compound system S + R, and the same argument following 
Eq. (27) holds. We therefore recover the Margolus-Levitin 
inequality for the purified compound system: 



(Kr) 



% h 

2jK) 



(43) 



The Mandelstam-Tamm and Margolus-Levitin inequalities do 
not hold for the most general quantum evolutions, described 
by completely positive, trace preserving maps (CP maps). To 
demonstrate this, we construct a counterexample. We can al- 
ways describe a non-unitary evolution of a single system as 
the unitary evolution of the system and its environment com- 
bined. The non-unitary evolution is obtained when we trace 
out the environment of the system. In our counterexample, 
the system starts and ends in a separable state with respect to 
the environment, which allows us to compare the states before 
and after the interaction. We show that orthogonal initial and 
final states can be created in arbitrary short times. 

Consider a two-level system (a qubit) with an energy gap Eq 
between the ground state |0) and the excited state When 
the qubit is prepared in the eigenstate of the Pauli X operator, 
the minimum time for this system to evolve unitarily to an 
orthogonal state is 



(45) 



saturating both bounds, with AE = {H) = Eo/2. In addition, 
consider qubits of a possibly different species with energy 
gap E^j. Our setup is shown in Fig. 1. We call the qubit with 
energy gap Eq the 'central qubit', and the remaining qubits 
are the 'satellite qubits'. At time t — 0, all qubits are prepared 
in the state |+) = (|0) + |l))/-\/2- To describe the counterex- 
ample we use the stabilizer formalism, in which the state of 
a system of qubits is fully determined as the +1 eigenvalue 
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FIG. 1 : Non-unitary evolution of the central qubit that violates both 
the Mandelstam-Tamm and Margolus-Levitin bound: (a) All qubits 
are prepared in the state |+); (b) each qubit is entangled simulta- 
neously with the central qubit using a CZ gate; (c) Hadamard gates 
are applied to the satellite qubits, which creates a GHZ state of size 
N+l; (d) phase shifts (j) are applied to the satellite qubits, and (e) 
a second set of Hadamard gates are applied to the satellite qubits; 
(f) A' simultaneous CZ gates disentangle the satellite qubits from the 
central qubit, leaving it with an accumulated phase shift of N(j> . 



N+l qubits. The phase shift on the satellite qubits given 
hy Uj{(j)) ~ exp(/0Z/2) change the stabilizer generators to 

N N 

4 ^^oYlUji<j))XjU]i<l)) =XoYl{cos<l)Xj + sm<j)Yj), 

(48) 

and the stabilizer generators for the satellite qubits remain un- 
changed: S^^^ = Sj^ . This evolution will take a time f^^^' . Af- 
ter the phase evolution, we again apply Hadamard gates to the 
satellite qubits, leading to the stabilizer generators 

N 

S'^^^ = Xo fj (cos Zj - sin (j) Yj ) and sf = ZoXj . (49) 

Finally, the satellite qubits are disentangled from the central 
qubit with another CZ gates, leading to 



r.(f) 



jQ -cos{N(j))Xo + sm{N<j))Yo and sf=Xj. (50) 

In other words, the satellite qubits are back to their initial state, 
and the central qubit is in the state 



|0)+e'^^|l) 

71 



(51) 



This is orthogonal to the initial state of the central qubit when 
N(j) = K. The total time T taken by this evolution is 



T = 2f, 



(Af) 
CZ 



-It 



H '^U ■ 



(52) 



In order to show that this is a counterexample to the 
Mandelstam-Tamm and Margolus-Levitin bounds, we need to 
show that we can choose parameters such that T < nh/Eo. 

The Hadamard gate evolves a state halfway to an orthog- 
onal state, and the speed limit for this gate can therefore be 
taken as half the time given in Eq. (45), with Eq ^ 



tH 



7th 

2E„ 



(53) 



state of a set of operators Si, called the stabilizer generators 
that generate an abelian group. For a detailed exposition of 
the stabilizer formalism, see Nielsen and Chuang [21]. The 
initial state of the central qubit (labelled '0') is then deter- 
mined fully by the operator 5g^^ =Xq. The satellite qubit j is 
stabilized by sf = X/, with ./ = 1 , . . . , A. 

Next, we entangle the qubits using CZ gates, causing the 
stabilizer generators to evolve to 



4^C=^on^7 =ZoX,-. (46) 



This evolution occurs in time f^'^' ■ The subsequent Hadamard 
gate on the satellite qubits changes the stabilizer generators to 



.7=0 

AN) 



(47) 



This evolution will take a time . It is straightforward to 
check that the state state corresponding to these stabilizer gen- 
erators is the GHZ state (|0,...,0) + |1, . . . , l))/\/2 on the 



Since all Hadamard gates are applied simultaneously to dif- 
ferent qubits, the minimum time to apply A Hadamard gates 
is the same as that for a single Hadamard gate. Therefore, we 
have f|^' =tf{- A similar argument holds for the application 
of the phase shifts on the satellite qubits, leading to r^^ — . 

The determination of the time it takes to apply the A CZ 
gates requires a little more care, since each gate involves the 
central qubit. First, consider the time evolution of Sj 



dt 



(54) 



The interaction Hamiltonian for a CZ gate between the central 
qubit and the satellite qubit can be written as Hj = gilo — 
Zq) (ly — Zj). Since — for all j and k, the interaction 

Hamiltonian for the A CZ gates is H = Solving the 

Heisenberg equation in Eq. (54) gives 



Sjit) = 



1 +COS 



1 — COS 



Xi 



1 



sin 



(!) 



ZqYj - - sin 



ZoYj,i55) 
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which exhibits a periodic behaviour with period T = 2nh/g. 
After half this period, the stabilizer Sj{t) becomes s''^\ and 
the gate time for a single CZ gate is therefore 



tcz • 



g 



(56) 



In the case of CZ gates, the evolution of the stabiUzer gen- 
erator So{t) is more complicated, since it encapsulates the in- 
teraction of the central qubit with all the satellite qubits. We 
have to show that it nevertheless exhibits the same periodicity 
as Sj{t). In general, we can write So{t) in terms of the Pauh 
operators on the central qubit and the j'^ satellite qubit, while 
collecting all other Pauh operators in Q^J^ : 



So{t) = ajit)Q'iiholj + Pj{t)Q'^'Yolj 

+7,(f)ei," ZoZ,- + 8jit)ZoQf YoZj , (57) 

where cij{t), jij{t), Yj{t), and 5j{t) are the coefficients that 
determine the time behaviour. Their exact form is not im- 
portant for our argument. Since the operators So{t) and Sj{t) 
stabiUze the evolving state at all times, and the stabilizer is an 
abelian group, we require that [S(){t),Sj{t)] = for all t. The 
periodicity of Sj{t) then imphes the periodicity of So{t), and 
moreover, they must have the same period. Since this period 
is independent of A^, we arrive at the conclusion that all CZ 
gates can be applied simultaneously. Notice that although the 
coefficients of Soit) wiU change for larger N, the underlying 
periodicity cannot be affected without causing the commuta- 
tion relations [So{t),Sj{t)] to become nonzero. 

Finally, we choose t(j, (and therefore (j) itself) arbitrarily 
small, and arbitrarily large, such that A^^ = n and J,^, w 0. 
The parameter regime where the system becomes a counterex- 
ample to the quantum speed Umits (t < nh/Eo) then becomes 



2nh 7th nh 

— + — < 



g 



(58) 



This can be re-cast into the following two requirements: 



Eg > Eo and g > 



■Eo 



(59) 



Since this is a perfecdy valid parameter regime in quantum 
mechanics, this constitutes a proper counterexample (even 



though it may be very difficult to implement in practice). It 
demonstrates that the bounds are valid only for unitary evo- 
lution. In fact, T can be made arbitrarily small by making Eq 
and g arbitrarily large, and consequently there is no general 
bound that is valid for all possible (unitary and non-unitary) 
quantum evolution. 

The violation of the Mandelstam-Tamm and Margolus- 
Levitin bounds leads to a number of open questions. First, 
what are the general requirements on a system and its interac- 
tion with the envirormient in order to violate the speed hmit? 
The counterexample presented here suggests that the interac- 
tion between the system and its environment must be strong, 
and that the environment must be more energetic than the sys- 
tem {Eg > Eo). Second, can these regimes be probed with 
experiments? It is not clear a priori that the strong inter- 
action of our counterexample can be achieved in a noiseless 
way. This may prevent the practical implementation of the 
evolution, even though the counterexample is perfectly valid 
from a mathematical point of view. In other words, taking 
into account noise may lead to another practical speed limit 
after all. And third, are there general bounds on the speed 
of quantum evolution based on general properties of the envi- 
ronment, rather than on detailed knowledge of the dynamics 
of the combined system? These questions must be addressed 
in order to fuUy understand the speed of dynamical evolution 
of open quantimi systems. 



V. CONCLUSIONS 

The Mandelstam-Tamm and Margolus-Levitin inequalities 
play an important role in the study of quantum mechanical 
processes in Nature. However, to date there has been only one 
derivation of the Margolus-Levitin inequality. In this paper, 
we gave alternative derivations for both inequalities from the 
statistical distance between quantum states. This allows for 
a fuUy geometrical interpretation of the quantum speed hmit. 
The inequalities were shown to hold for unitary evolution of 
pure and mixed states, and a counterexample to the inequal- 
ities is given for evolution described by completely positive 
trace-preserving maps. 

The authors thank Carlos Perez-Delgado, David Whittaker, 
and Samuel Braunstein for valuable discussions. 
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